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The spinor Bose-Einstein condensate of atomic gases has been experimentally realized by 
a number of groups. Further, theoretical proposals of the possible vortex states have been 
sugessted. This paper studies the effects of the quadratic Zeeman energy on the vortex states. 
This energy was ignored in previous theoretical studies, although it exists in experimental 
systems. We present phase diagrams of various vortex states taking into account the quadratic 
Zeeman energy. The vortex states are calculated by the Gross-Pitaevskii equations. Several 
new kinds of vortex states are found, ft is also found that the quadratic Zeeman energy affects 
the direction of total magnetization and causes a significant change in the phase diagrams. 

KEYWORDS: Vortex, Spinor, Bose-Einstein condensate 



1. Introduction 

Quantized vortices in the Bose-Einstein condensate 
(BEC) of ultracold atoms have been a topic of active 
research. Since the first realization of the vortex, 1 a vari- 
ety of studies have been carried out such as those on the 
formation of the vortex lattice, bending and oscillation 
of vortex lines, 2 vortex rings, and so on. For a review of 
the early studies, see ref. 3. The interest in vortices in 
ultracold atoms is also increasing with regard to more 
complicated systems, such as the BEC interacting with 
an optical lattice, BEC with dipolar interactions, 4-6 and 
molecules and Cooper pairs of Fermions. 7 

Despite these extensive studies on vortices, studies on 
simpler systems such as the spinor BEC of atoms in a 
uniform magnetic field 8-10 remains unachieved. The ex- 
istence of various vortices has been predicted in the sem- 
inal papers by Ho 11 and Ohmi and Machida 12 and was 
subsequently continued by many works. 8 ~ 10, 13-15 

In the spinor BEC, a vortex state is inseparably related 
to the spin texture. By utilizing this relation, it is pos- 
sible to imprint vorticity without rotating the cloud as 
proved in the experiments. 16 By applying a quadrupole 
(nonuniform) magnetic field, the spin texture of the BEC 
is forced, resulting in the nucleation of the vortex. 

Conversely, for a rotated spinor BEC in a uniform 
magnetic field or without a magnetic field, many the- 
oretical predictions on the possible spin textures have 
been made. 8,9 ' 13-15 In this case, the rotation of the trap- 
ping potential causes vorticity that produces a spin tex- 
ture. In a spinor BEC in a uniform field, the boundary 
condition and magnetic field gradient do not explicitly 
control the vorticity. Instead, the internal energy of the 
vortices determines the vortex configuration. However, 
the experimental observation of vortices still remains to 
be performed. 

In an optically confined BEC, a dominant effect - 
the linear Zeeman effect — due to the intensity of the 
magnetic field is cancelled by the conservation of the to- 
tal magnetization (see ref. 17). Nevertheless, the mag- 



netic field intensity affects the internal spin through the 
quadratic Zeeman effect, which was ignored in previ- 
ous studies 8,9 on the vortices in the spinor BEC. The 
quadratic Zeeman energy is proportional to the square 
of the magnetic field intensity and the local spin. 

In this paper, we take into account the quadratic Zee- 
man energy and present new phase diagrams of vortices 
calculated by the extended Gross-Pitaevskii equation. It 
is important to consider the quadratic Zeeman energy 
because its existence in reality is inevitable. In addition, 
the quadratic Zeeman energy leads to vortex reconfigura- 
tion within fixed values of magnetization M and angular 
momentum L. 

Atomic interactions are classified into two types: fer- 
romagnetic and antifcrromagnetic (polar). 11, 12 In this 
study, we treat the condensate of sodium atoms (F = 1), 
which is antiferromagnetic. In the absence of a magnetic 
field for antiferromagnetic systems, vortices (and their 
lattices) in which the up and down components are pop- 
ulated are stable. 13,15 They have been labelled as the 
vortex of region II 8 or the (1, x,0) type. 9,10 In a mag- 
netic field, or in the presence of linear Zeeman energy, 
various other vortex states are found to be possible. The 
present work aims to study the effects of the quadratic 
Zeeman energy on such vortex states and to extend the 
phase diagrams. 

2. Model 

We simulate the vortex states of the BEC of atomic 
gases by numerically solving the extended Gross- 
Pitaevskii equation. The system is rotated with an an- 
gular velocity f2, and a uniform magnetic field is applied. 
The calculation is similar to that in ref. 8, but we take 
into account the quadratic Zeeman energy. For simplicity, 
we treat a harmonically trapped two-dimensional system. 

We assume the direction of the applied magnetic field 
and the direction of quantization axis to be the z-axis. 
The vector of the rotation f2 is also in the +z direction. 

The order parameter of the condensate consists of 
macroscopic wavefunctions where j — 0, ±1 are hy- 
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pcrfine projections. The order parameters iffj are given 
by the minimization of a free energy 

G = E 1 + EU + E Zcc — flL — nN, (1) 

which consists of the single particle energy E\ , the inter- 
action energy E mt , and Zeeman energy E Zee ; they are 

E 1 = ]T / dr*J (-^V 2 ) ^ + V(r) nj (2) 

E Zcc = ^2 dr (-pjrij + qfrij) , (4) 

respectively, where the densities rij = \^fj\ 2 and n = 
rij, particle number N = . J" rij(r)dr, angular mo- 
mentum i = —ih^2jf ^*V^>j x rdr, spin index j e 

1,0,-1, and trapping potential V(r) = m ( 2 ^) 2 (x 2 +y 2 ). 
The term proportional to p (q) in eq. (4) represents the 
linear (quadratic) Zeeman energy. 

The minimization of G gives the established Gross- 
Pitaevskii (GP) equation for the spinor BEC. 11,12 We 
solve the GP equation numerically for various values of 
the quadratic Zeeman factor q to investigate its effect on 
the vortex configurations. The parameters p and p are 
varied to obtain results with various magnetization M z 
at fixed N = N 2 d- 

The particle number of the j-th component is Nj = 
J rij(r)dr. The normalized angular momentums of the 
j-th component Lj = — ifiJ2j I *j x rdr/Nj where 
the spin index j e 1, 0, —1 arc extensively used to identify 
vortices. The local and total magnetizations are m a (r) = 
\I/t(r).F a \I>(r) and M a = J m a (r)dr, respectively, where 
a = x,y,z. 

A system parameter e is defined as follows. From the 
mass M a of atom, an effective scattering length a = 
(do + 2a 2 )/3, the interaction parameter Co = Aith 2 a/M a 
and harmonic trap frequency u>o, we can calculate a 
harmonic oscillator length Ao = y fyMaMJo- We use 
a ratio of oscillator length to condensate radius -Rtf, 
£ = (Ao/-Rtf) 2 as a system parameter. Here the ra- 
dius Rtf is that of condensate within Thomas-Fermi 
(TF) approximation in the absence of spin interaction. 
For the system parameter, values e = 0.1 and 0.01 are 
used. After a e is given, we have a particle number 
^2D = l/(16e 2 a) within the TF approximation. We 
fix the particle number N to A^d- The approximate 
peak density becomes n pea k = hu r /(2coe) where h is 
Plank's constant. For sodium atom when v r = 100 Hz, 
n = 0.329 (3.29) x 10 14 cm~ 3 for e = 0.1 (0.01). 

Value e = 0.1 is used for consistency with our previ- 
ous study. 8 Another value e = 0.01 is employed to sim- 
ulate thin disk-shaped condensate. If we assume tight 
harmonic confinement and Gaussian shape of conden- 
sate wavefunction along z-axis, the particle number in a 

system is A 3 d = \J ^2D where uj z is the confine- 
ment parameter along z-axis. For sodium atoms when 



uj z = 1000 x 2tt and e = 0.01, iV 3D is 2.65 x 10 5 . This is 
within typical experimental particle number. 

In an experiment 18 without a vortex, q ranges from 
to 70 Hz for magnetic fields B up to 500 mG. We 
choose q/{hv r ) from to 0.1. The maximum 0.1 cor- 
responds to 190 mG when v r = 100Hz and q oc Bq. 

The parameters that we explicitly give its values are 
the system parameter e, the angular velocity O and the 
quadratic Zeeman factor q. The total magnetization M z 
and the particle number TV are controlled via the linear 
Zeeman factor p and chemical potential fx. We employ the 
ratio of the interaction parameters C2/C0 = 0.02, which 
is obtained from an experiment 18 that directly treats 
the quadratic Zeeman energy and interaction parameter. 
Similar to that experiment, we assume the condensate of 
sodium atoms with F — 1. 

There are a number of vortex states at each M z . 
Among them, we assume that a state with the lowest 
energy G is realized. These lowest energy states are pre- 
sented in the phase diagrams — Figs. 1(a), 1(b) and 4. 
However, the direction of the total magnetization re- 
quires further explanation before the details of the phase 
diagrams are examined. 

3. Phase diagrams 

3. 1 Direction of total magnetization 

Total magnetization is conserved in an optically con- 
fined condensate (see refs. 17, 18). For this conserva- 
tion, the factor p of the linear Zeeman energy —pM z , 
which was originally proportional to the magnetic field, 
is treated as a Lagrange multiplier. Therefore, the fac- 
tor p loses its direct correspondence with the magnetic 
field. Conversely, the quadratic Zeeman factor q, which 
is proportional to the square of the magnetic field, is not 
restrained by the conservation. 

The horizontal axis in the diagrams is the normalized 
magnetization M z /N in the z-direction. The total mag- 
netization M x , My, M z may be rotated around z-axis by 
the spin rotation exp~ lF * r , without affecting the Zeeman 
energies. Therefore, elements M x and M y are denoted by 

Mi = ^Ml + M 2 hereinafter. 

The amplitudes of the elements m x and m y of the local 
magnetization, and therefore the amplitude of Mj_, arc 
proportional to the amplitude |^o| of the zero species. 
Because the zero species of the condensate becomes en- 
ergetically favorable for larger q, the total perpendicular 
magnetization M± can be large for some types of vor- 
tices. 

Usually, the total magnetization M± perpendicular to 
the z-axis is zero. However, in some regions, vortices with 
M± have the lowest energy. The regions containing 
vortices with M± ^ are also shown in Figs. 1(b) and 4 
as the area shaded with (blue) slant lines. 

3.2 Absence of quadratic Zeeman term 

In the absence of the quadratic Zeeman energy (q = 0), 
our system is equivalent to that in ref. 8 at C2/C0 = 0.02. 
In Fig. 1(a), regions IV, II and III, which have been re- 
ported before, 8 are observed. The boundaries of these 




Fig. 1. Phase diagrams of magnetization M z /N vs the quadratic 
Zeeman term parameter q/(hu r ) at e = 0.1 and f2 = 0.45. (a) 
This diagram shows the areas for the vortices with M± = 0. A 
set of three numbers in parentheses e.g., (1, 0, —1) represents the 
angular moments L\ , Lo and L_i of the condensate, respectively. 
In a shaded region, when the numbers in the set are integers, all 
the vortex states in that region has that set of angular moments. 
Otherwise, a set of numbers shows only one typical value in the 
region. Letters u, 0, or d in double circles indicates that the 
up, zero, or down species is dominant inside (inner circle) or 
outside (outer circle) the centered vortex core, (b) Regions II-l 
and II-2 in which vortex states with M± ^ have a lower energy 
are shown. The perpendicular magnetization M±/N in II-l lies 
between 0.001 ~ 0.2 at qjv r > 0. At q = and M z /N > 0.5 
within II-l, Mx/N is as low as 10~ 8 . 

regions are consistent with the present results. However, 
there are also vortex regions II-l and II-2 exist with 
Mj_ 7^ and a region III-l also exists with M± = 0. 
Figures 1(a) and 1(b) show phase diagrams along the 
horizontal axis (q = 0). Let us consider these regions 
from smaller M z values at Mj_ = 0. Amplitudes of 



of the corresponding vortices are shown in Fig. 2. 

IV, IV- 1: The zero species has a vortex (Lq ~ 1) and 
the other species fill the core region. Moreover vortex 
IV- 1 has integer angular momentums (Li,Lq,L-{) = 
(0,1,2). For a simple explanation, let us consider the 
cases with a single centered vortex on one or two of the 
three components, and that the other component (s) fillx 
the vortex core. For a smaller magnetization M z < 0.2, 
the local magnetization m z /n = (ni— n_i)/n tends to be 
zero. When the local magnetization m z /n is zero, there 
are two possible configurations: (1) There is a vortex 
in the zero component, and a mixture of up and down 
species fills the vortex core, and (2) the mixture has a 
vortex (up and down species are co-rotating), and the 
zero species fills the core. Here, under the phase condi- 
tions introduced below in eq. (7), (1) is preferred. Hence, 
type IV is realized. As M z is increased from zero, the an- 
gular momentum L—± smoothly increases to fully satisfy 
eq. (7) in the IV- 1 area. 

III-l: As mentioned above (2), the up and down species 
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exhibit vorticity L\ and L_i <~ 1, and the zero species 
fill the core. 

II and III: The total angular momentum around L = 
0.7 is preferred due to the angular velocity ft. Therefore, 
the dominant species should have one vortex. For a larger 
M z (> 0.2), throughout the II and III regions, the up 
component is dominant and has one vorticity L\ = 1. In 
region II (III), the down (zero) component fills the core 
to decrease (increase) M z . The type II and III vortices 
correspond to (1, x,0) and (1,0,-1) in rcf. 10, respec- 
tively. In this notation, the (integer) angular moments of 
each component are listed in parentheses for the axisym- 
metric cases. 

Let us consider the vortices with finite M± as shown 
in Fig. 1(b). As compared to the corresponding vortex 
states in Fig. 1(a), the vortex states in Fig. 1(b) have a 
lower energy G. The II- 1 and II-2 regions are shown in 
Fig. 1(b), and these two types of vortices can be smoothly 
deformed into each other. 

From M z /N = to 0.276: For specific value of the 
parameter p — 0, a configuration of two vortices is pos- 
sible. This includes the "split- (I)" states in ref. 10. In 
our case, the vortex state has M z /N — and M±/N = 
0.276. Because the Zeeman factors are zero, the direc- 
tion of spin may be freely rotated by the operators 
exp -4 ^" cxp~ lFy P. By this rotation, this vortex is re- 
duced to a type II vortex with M±/N = 0,M z /N = 
0.276. By similarly continuing the rotation, these vor- 
tices cover a range of M z /N = — 0.276 while main- 
taining \/M\ + Ml /N = 0.276. Similar deformations 
of the density profiles among the spin components due 
to changes in the quantization axis are reported in 
refs. 8, 15. This deformation affects the energy when the 
Zeeman term factors p, q are finite, and thus no special 
care is required in other areas of the phase diagrams. 

From M z = 0.51 to 0.66: The vortex states in this 
range have a small finite perpendicular magnetization 
M±/N <~ 10~ 8 , and these are similar to the vortices 
within < M z < 0.276. Nevertheless, these cannot be 
reduced to vortices with M± = because of the finite 
Zeeman term factor p. This region is discussed in §3.4. 

3. 3 Effects of quadratic Zeeman term 

Let us consider the effects of the quadratic Zeeman 
energy shown in eq. (4). This energy is proportional to 
the parameter q. An increase in the Zeeman energy for 
a finite q is q{N\ + -/V_i) = q(N — No). Therefore, there 
is a maximum increase in the energy of vortex of type 
II, whose N — 0; it becomes less favorable. In the phase 
diagram, the type II region shrinks above q/(hu r ) = 0.05, 
and a large range of M z is occupied by the other regions. 
The other boundaries of the regions are not sensitive to 

The non-axisymmetric region III-l exists at a finite q. 
In region III-l, angular moments L\ and L_i have the 
same sign and a similar amplitude. In other words, the 
up and down species are co-rotating in region III-l, while 
they are counter-rotating in region III. The c 2 term of 
the energy prefers the counter-rotating configuration, see 
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Fig. 3. (a)(b) A vortex state within the region II-l at e = 0.01, 
q = 0.06hu r , M-JN = 0.437 and M ± /N = 0.080. (c)(d) A vortex 
state within the region III-2 at e = 0.01, q = 0.05hi/ r , M z /N = 
0.302 and M±/N = 0.085. (a)(c) Spin textures. Arrows show 
the directions of local spin components m x and m y . Component 
nriz/n is shown by contours. (b)(d) Phases of V&o- Solid lines 
indicate multiples of 7r/4. 



discussion in §3.4. However, the co-rotating configuration 
becomes energetically favorable by decreasing the energy 
—uuL because the total angular momentum L <~ {N\ + 
N-^/N is larger than L = (Ni—N_i)/N in the counter- 
rotating case. 

Several data points show the III-l configurations even 
within the region II on Fig. 1 (a) . Their energies are lower 
than the type II vortices. But we couldn't make sure 
that they (III-l) always exist within certain range of M z 
because they become unstable by small changes in p. 

3.4 Phase in the region II- 1 

In the region II- 1, the phases and angular momen- 
tums of the up and down species are similar to those in 
the region II. However, the wavefunction is nonzero 
there. This minor change in the wavefunction signifi- 
cantly changes the spin texture because the components 
of local magnetization m x and m y are proportional to 
the amplitude |^o|- The spin texture in the xy plane ap- 
pears as shown in Fig. 3(a). The phase of <&o is important 
for understanding of this spin texture. 

The phase of ^0 is restricted by the minimization of 
the local spin interaction energy of Ei nt in eq. (3). The 
C2-proportional part of the local energy is written as 

S c2 = In+ltn 2 -^!*-!-^!) (5) 

= C + 2c 2% /nin_in cos(2a — ct\ — a_i) (6) 

where ctj are the phases of "J/j, and C is a function of 
densities nj . In our antiferromagnetic system (c 2 > 0) , a 



condition for minimizing £ c2 for given rij is 

2ao — ol\ — a_i = 7r + 2m7r, 
where m is an integer. 



(7) 
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Fig. 4. (Color online) Phase diagram of magnetization M/N vs 
quadratic Zeeman term q/(hv r ) at e = 0.01 and Q = 0.07. No- 
tations are same as in Figs. 1(a) and 1(b). But both Mj_ = 
regions and M± ^ regions (shaded with slash lines) are shown 
here. A region at g = 0.07, M z /N ~ 0.25 is not labelled. There 
angular momentums are around (2.7,0.3,0.3). 



In regions II and II- 1, a\ ~ 9, a_i = 0. Therefore, 
ao = | and the zero species must have a half wind- 
ing number. In region II-l, instead of a half winding 
number, the zero species has an off-centered vortex core 
and a non-integer angular momentum < Lq < 1, see 
Fig. 3(b). 

The components of the magnetization m x and m y are 
proportional to the amplitude \^o\. In II- 1, because the 
amplitude \^q\ becomes finite, the spin texture in the xy 
plane appears as shown in Fig. 3(a). 

In the absence of Zeeman terms p = q = 0, this struc- 
ture may be reduced to vortices in region IV by changing 
the quantization axis. At a finite p or q, these vortices 
cannot be reduced similarly and are energetically more 
favorable than type II vortices. Therefore, the II-l vortex 
is likely to be observed in experiments conducted with fi- 
nite magnetic fields. 

3.5 Disk-shaped condensate 

Figure 4 is the phase diagram as e = 0.01. Amplitudes 
of YS>i\ of the corresponding vortices are shown in Fig. 5. 
This value corresponds to the realistic particle number 
for tight z-confinemcnt traps. We employ the angular 
velocity £1 = 0.07 x 27r to cover the vortex states observed 
at e = 0.1. At a lower (higher) rotation, = 0.05 x 2tt 
(0.1 x 2tt), vortexfree (many vortex) states replace these 
vortex states. 

At q = 0, regions I, II and III are observed. Region I 
is reported in ref. 8. The composite vortex in this region 
consists of two vortices, one is in the up component and 
the other is in the down component. The zero compo- 
nent does not have an atom. The composite vortex in 
region I at p = q = is equivalent to that in region IV 
through the tt/2 rotation of the quantization axis. The 
other regions III and IV are similar to the abovemen- 
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Fig. 5. Amplitudes of (i = 0,±1) of vortices in regions I, II, 
III, III-l, IV and IV-1 at e = 0.01 and Q = 0.07. The normalized 
angular momentum Li is shown in each panel. 

tioncd e = 0.1 case. 

At finite q, because the absence of the zero component 
is energetically unfavorable, region I is replaced by region 
III-l. The vortices in region III-l may be understood as 
those in region I with finite atoms of zero species. 

Regions II, III, III-l and IV are similar to the e — 
0.1 case. However, the region II elongates toward larger 
q because the scale of e changes with q. Consider two 
energy densities at the trap center E c i ~ C2n 2 oc C2£~ 2 
and E q ~ qn oc <7£ _1 , where n is the atom density. These 
energies determine the relative population of the zero 
species, and they should be comparable at the boundary 
of region II. When we consider the boundary between 
systems with different e values, systems with equal q£ 
rather than q are in correspondence. 

In regions III, (2,0,-2) and (3,0,-3), the up and 
down species are counter-rotating with respect to each 
other (angular moments L\ and L_i have opposite sign 
but the same amplitude). These satisfy the phase rela- 
tion given by eq. (7). In region III-l, this condition is not 
met. Region III-l, in which the up and down species are 
co-rotating, lies in between regions IV and II. 

Vortices with finite perpendicular magnetization M± 
reduce the energy at larger p. They are shown as areas 
shaded with oblique lines, e.g., II-l and III-2 regions. 
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The II- 1 region at M z /N ~ 0.5 that also appeared when 
e = 0.1. Here, M±/N <~ 0.03 - 0.09. Another region is 
III-2, in which M±/N ~ 0.07. A pair of spin structures 
similar to those of the type II- 1 vortices exist at e = 0.1. 

When the initial perpendicular magnetization is zero 
(M±/N = 0) and the total magnetization is explicitly 
conserved, only the vortex configurations with M±/N — 
are selected. If there is an error in M±/N, the II- 1 or 
III-2 type of vortices having lower energy will be real- 
ized. Therefore, determining whether the areas shaded 
by oblique lines should be considered in this phase dia- 
gram depends on the experimental details. 

The vortices in region III-2 have a pair of the line- 
like phase singularities on the zero species as shown in 
Fig. 3(d). Figure 3(c) presents resulting spin texture. Be- 
tween the singularities, the center of the trap is filled by 
zero species. These vortices are regarded as a pair of type 
II vortices that has gained finite zero species ^o- 

4. Conclusion 

We investigated the effect of the quadratic Zeeman 
energy on various vortex states of an antiferromagnetic 
spinor BEC. We presented phase diagrams showing the 
quadratic Zeeman term factor q vs magnetization M z /N. 
An important observation is that the perpendicular mag- 
netization M±/N is also necessary to correctly plot the 
phase diagram, not only at finite q but also at q = 0. The 
permissible range of M± will depend on the experimental 
conditions. 

In a specific case (p = q = 0) in which the Zeeman 
terms are ignored, the rotation of all the spins deforms 
a vortex in region II into other configurations (regions 
II-l and II-2) with a finite perpendicular magnetization 
Mj_. This cov ers the ran ge < M z /N < 0.276 while 
maintaining y/M% + M\/N = 0.276. 

For a finite q, due to the quadratic Zeeman energy, 
vortices with M± ^ (in regions II- 1, II-2 and III-2) 
widely replace the other vortices and change the phase 
diagram. In the general case (p 7^ or q 7^ 0), configura- 
tions with finite M± cannot be reduced into configura- 
tions with M± = 0. Axisymmctric vortices (3, 0, —3) and 
(2,0,-2) whose up and down components are counter- 
rotating with respect to each other are also found at 
larger q. 

In the absence of the quadratic Zeeman term (q = 



0), the result is consistent with previous studies by one 



of 



,c 8 



However, vortices II-l and III-l are also found. 



These vortices may be understood as a deformation of 
the vortices in regions II and I with finite zero species. 

The vortices in region II are commonly known for 
the antiferromagnetic spinor BEC. 10 ' 13 ' 15 For a finite q, 
these vortices gains zero species and are deformed into 
II-l vortices. Similarly, III-2 vortices are regarded as de- 
rived from a pair of type II vortices. A similar deforma- 
tion of vortices by having finite zero component is ex- 
pected in the spinor BEC in a magnetic field, especially 
in the lattices of the type II vortices. 15 
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